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FOREWORD 


The  fi lamentation  instability  properties  of  a  relativistic  hollow  electron 
beam  confined  in  axial  flow  by  a  uniform  magnetic  field  in  a  pipe  are  investigated 
via  the  VI asov-Maxwel 1  equations.  The  instability  is  found  to  have  two  sidebands, 
one  with  a  spectrum  of  positive  wavenumbers  k  and  the  other  with  a  spectrum  of 
negative  wavenumbers.  The  spectral  point  k=0,  associated  with  the  diokotron 
instability,  is  excluded  from  the  fi lamentation  instability's  two  unstable 
sidebands.  Only  in  the  limit  of  zero  axial  beam  flow  (y-*-l),  the  diokotron 
instability  becomes  asymptotically  part  of  the  filamentation  instability  spectrum. 
In  this  limit  the  filamentation  instability's  two  sidebands  merge  asymptotically 
and  symmetrically  toward  the  diokotron  instability  spectral  point,  k=0,  in 
agreement  with  the  basic  driving  physical  mechanisms  and  geometry  configurations 
for  these  two  distinct  and  different  instabilities. 

IRA  M.  BLATSTEIN 
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INTRODUCTION 


In  a  recent  paper  by  Uhm  and  Siambis*  the  diokotron  instability  of  a 

hollow  relativistic  electron  beam  in  a  conducting  pipe  guided  by  a  uniform  axial 

magnetic  field  was  investigated.  Relativistic  and  electromagnetic  effects  were 

included  in  the  derivation  of  the  properties  of  the  instability.  In  Reference  1 

the  analysis  of  the  diokotron  modes  proceeded  by  taking  the  limit  k=0,  ab  initio, 

where  k  is  the  wave  number  along  the  beam.  This  restriction,  of  k=0,  for  the 

diokotron  modes  has  been  assumed  in  all  earlier  treatments  of  the  diokotron 
2-4 

instability  .  Physically  this  assumption  was  motivated  by  the  geometrical 

2 

configuration  of  the  magnetron  tube  in  particular,  as  well  as  other  crossed 

3  5 

field  (E-cross-B)  beam  devices  ’  utilizing  hollow  electron  beams.  Instabilities 
in  hollow  beams,  in  long  cylindrical  pipes,  confined  in  axial  flow  by  a  uniform 
magnetic  field,  were  investigated  experimentally  by  Kyhl  and  Webster®  and 
analytically  by  Pierce7  in  order  to  establish  the  mechanisms  for  noise  generation 
and  amplification  in  traveling  wave  microwave  tubes  with  hollow  beams.  These 
early  investigators6,7  found  that  low  voltage,  low  current  hollow  electron  beams 
can  break  up  into  vortex  filaments,  similar  to  those  associated  with  the  diokotron 
instability,  and  in  addition  they  exhibit  a  non-zero  wave  number,  k/0,  along  the 
applied  magnetic  field.  They  observed  that  the  azimuthal  wave  length  of  the 
instability  was  comparable  to  the  longitudinal  wave  length  of  the  instability. 

They  also  observed  that  this  comparabil ity  in  axial  and  azimuthal  wave  lengths 
was  maintained  as  the  instability  actually  evolved  from  the  fastest  growing  short 
wave  lengths  to  the  slower  growing  longer  wave  lengths.  They  called  this  in¬ 
stability  the  filamentation  instability.  These  early  investigators  were  interested 
in  whether  this  phenomenon,  at  low  levels,  contributed  to  microwave  tube  beam 
noise.  They  concluded  that  a  direct  contribution  to  beam  noise  seemed  unlikely 
since  the  waves  did  not  have  the  proper  symmetry  to  couple  to  the  rf  circuits 
of  the  microwave  tubes®.  This  observation  resulted  in  lack  of  broad  further 
interest  in  this  phenomenon. 
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More  recently  intense  relativistic  hollow  electron  beams  have  become 
the  object  of  intense  experimental  and  analytical  investigations  in  connection 
with  a  broad  spectrum  of  modern  applications  such  as  high  current  electron  beam 
accelerators,  collective  accelerators,  gyrotrons,  free  electron  lasers  and  fusion. 
It  is  the  purpose  of  this  work  to  investigate  the  properties  of  the  filamentation 
instability  for  intense,  relativistic  hollow  electron  beams  of  interest  in  current 
applications.  In  carrying  out  the  analysis  we  shall  follow  the  technique  of 
Reference  8.  In  Section  II  the  relativistic  electron  hollow  beam  equilibrium 
state  is  analysed.  In  Section  III  a  Vlasov-Maxwell  stability  analysis  is  carried 
out.  In  Section  IV  the  axisymmetric,  1=0,  space  charge  mode  for  the  hollow  beam 
is  obtained  and  found  to  be  in  agreement  with  results  from  more  approximate 
theories.  In  Section  V  the  filamentation  instability  modes,  1>3,  are  analysed 
and  discussed.  In  Section  VI  conclusions  are  presented. 
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equilibrium 


We  consider  the  equilibrium  configuration  illustrated  in  Figure  1, 
consisting  of  an  intense  hollow  relativistic  electron  beam  propagating  in  a 
drift  tube  parallel  to  a  uniform  applied  magnetic  field,  BQ§Z>  with  velocity 
^bc^z‘  Cylindrical  coordinates  (r,9,z)  are  used,  with  the  z  axis  along  the  axis 
of  symmetry.  The  beam  is  described  by  a  distribution  f(r,p,t)  which  satisfies 
the  Vlasov  equation 


H 


f(r,p,t)  =  0 


(1) 


where  r  stands  for  the  cylindrical  coordinates,  the  momentum  p  for  p  ,pa,p  , 

I  w  A 

the  fields  £  and  §  are  the  external  and  self-consistent  fields  associated  with 
the  beam  flow.  The  velocity  v  is  given  by 


v  =  p/ym  ,  8  =  v/c 


(2) 


where 


and  -e,  m  are  the  charge  and  rest  mass  of  the  electron,  and  c  is  the  velocity 
of  light. 


The  hollow  beam  equilibrium  flow  distribution,  suitable  for  the 
application  at  hand,  is  given  in  terms  of  the  constants  of  the  motion. 

fb '  ^5  -  ebcpz  -  j(p 

where  the  total  energy 


V 


(4) 


H  =  (m2c4 


+  p2c4 ) 2  -ebQ(r) 


(5) 
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the  canonical  angular  momentum 

pe  ’  r<»9  -  H  Bor)  • 
and  the  axial  cononical  momentum 


are  the  three  single  particle  constants  of  the  motion. 

The  constants  nQ,  Rq,  Pq,  Yb  and  y0  are  identified  as  follows: 
quantity  nQ  is  the  value  of  the  density  nb(r)  at  the  median  radius  RQ  as 
in  Figure  2,  PQ  is  the  value  of  the  canonical  angular  momentum  at  r=RQ, 
is  the  total  electron  energy  at  the  beam  frame,  and  Yb  is  given  by 

yb  Ml  -  • 

We  also  define  the  energy  variable  U 
U  ■  H  -  BbcPz  -  YQmc2/Yb 
and  write  the  equilibrium  distribution  as 

fb  -  5<u>4<pe  *  V  • 

After  a  straightforward  algebra,  U  in  Equation  (9)  is  expressed  as 

U  =  mc2  -  e  [<t>Q(r )  -  SbAz(r)]J 

where  y'nc2  is  the  electron  kinetic  energy  at  the  beam  frame. 


(6) 


(7) 


The 

shown 

omc2 

(8) 


(9) 


(10) 


(ID 
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The  equilibrium  electromagnetic  fields  are  obtained  from  the  scalar 
potential  $  (r)  and  vector  potential  A  (r)  from  Maxwell's  equations 


?l?rlF4o('')  *  4,enb(r) 


?  If  r  If  \(r>  3  *Wr) 


02) 

(13! 


where  nh(r)  =  /  d3pfh(H,P0,Pz) 


(14) 


which  is  found  to  be 


nb(r)  = 


0 

nc 

0 


noRo/r 


r  <  Ri 
Ri  <  r  <  R2 
r  >  R2 


(IB) 


Equations  (12)  -  (15)  give 


Er(r) 


Be(r) 


34>„(r) 

-1F~-  -4^0Ro  (1  *  R</r> 
3A  (r ) 


R,  <  r  <  R: 


3r  =  -4ireBbnQRo  (1  -  Ri/r)  ,  Rj  <  r  <  R2 


(16) 

(17) 


The  radial  force  balance  on  the  beam  electron  field  element,  obtained  from 
Equations  (1)  and  (4),  is 


VV0(r) 


Er(r)  + 


VQ(r) 


Bo  '  6bB9(r) 


(18) 


We  define  the  equilibrium  rotation  frequency 


we(r)  =  VQ(r)/r 


(19) 


n 


V-  -  «t  iSTV  3*-\  ‘  V  2 1  «•  .  r 
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and  from  Equation  (18)  write  down  its  value  for  the  median  electron  fluid  element 


at  r=R 


o’ 


where 


eBo 

M  =  -  , 

c  mcYb 
47te2n, 

2  _  D 

(jJ.  =  - 

pb 


(21) 

(22) 


and  a  =  R2  -  R  =  R  -  R ,  is  the  half  thickness  of  the  hollow  beam.  The  equilibrium 
oo 

state  of  interest  is  the  one  with  the  slow  rotational  frequency  u)~(Rq).  Additional 
properties  of  the  general  equilibrium  state,  assumed  here  through  Equation  (11), 
are  similar  to  those  derived  in  Reference  8  for  the  non-relativistic  case. 
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STABILITY  ANALYSIS 

In  this  section  we  use  the  linearized  VI  a  so  v -Max'.-, 'el  1  equations  to 
investigate  the  general  normal  modes  (t  >  0,  k  f  0)  and  their  stability  pro¬ 
perties  for  the  type  of  intense,  hollow  beams  discussed  in  the  equilibrium 
section.  We  adopt  a  normal  mode  approach  in  which  all  perturbed  quantities  are 
assumed  to  vary  with  r,  9,  z,  and  t  as 

6vp( r , t)  =  <5^(r)  exp  [  i  (*.  8  +  kz  -  wt)j  .  (23) 

For  the  electromagnetic  fields  we  assume  the  TM  (transverse  magnetic)  modes  which 


have  the  following  field  components 

&t  =  +e  6E  +  efi6Ea  4  e  5E_  (24) 

r  r  6  9  z  z 

6§  *  §r5Br  4  e95B0  (25) 

and  which  satisfy  the  linearized  Maxwell's  equations 

$-6t  =  4tt6p  (26) 

Vx6?  =  ^  (27) 

Vx6B  =  —  <SJ  -  —  6t  (28) 

c  c 


Sbustitution  of  Equation  (28)  into  Equation  (27)  and  rearrangement  with 
Equation  (26)  generates  the  wave  equation  for  6E2 

5EZ  =  4irik(6p  -  oJ2)  .  (29) 

However,  the  axial  component  of  the  perturbed  current  density  is  approximated  by 


a_ 

3r 


(k2  -  UJ2 / c 2 ) 
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SJ2  =  8bc6p  .  (30) 

The  analysis  will  now  be  restricted  to  the  cases  of  interest,  where 
VQ/c  «  1  and 


With  these  assumptions  the  5Jr  and  <$JQ  source  terms  in  Equation  (30)  can  be  neglected. 
Also  for  wave  numbers  k  satisfying 

(k2R*/Y2)  «  1  (35) 

the  third  term  in  the  left  hand  side  of  Equation  (29)  can  be  neglected,  resulting 
in  the  approximate  wave  equation 


Next,  the  source  term  <$p  will  be  evaluated  by  first  finding  the  perturba¬ 
tion  <5f,  of  the  distribution  function  f^,  from  the  linearized  Vlasov  equation 
t 

6f  =  ef  d t ‘  exp(-iwt' )  exp(ikz  +  i £0)  6t  +  fb  (37) 

- 00  3p' 

where  the  integration  in  t'  is  carried  along  the  equilibrium  particle  trajectories. 
We  integrate  by  parts  with  respect  to  t'  and  change  variables  to  t  =  t'  -  t,  to 
find 
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(  9fh 

3fb 

5fh 

}3£2(r>  3U  *  1 

^  "  k£bc^  W  + 

0 

0 

where  the  orbit  integral  I  is  defined  by 
/•0 


I  * y* dT6Ez(r')  exp  -iurr  4  iJt(e*  -  e)  +  ik(z‘  -  z) 

—  oo 

In  order  to  evaluate  Equations  (37)  -  (39),  we  use 


e'  s  ur  „  -  u 


<P6  -  Po> 


eo  “  y.mRI 


u  =  u£/u>*  -  1 


2  .  [  2Po  \  ,  / .  4uipba  \ 

%  -UbmR0j  '  "cl1  ‘  ffe)  ■ 

We  assume  that  i  (u  /u  )  «  R  /a  and  we  approximate  6E  by 

W  J)  U 


P.  -  ?n  u>  /  3<$E  \ 

KZW)  -  -|(V)R 


so  that 

I 

where 


Pa-Prt  Un  /36E  \ 

=  1[6Ez(Ro)  + 


pz  *  yJ;  (PZ  -  Ybm0bc) 


“*  -  kBcc  '  1  (weo  '  u 


pa-p„  \  kp;  * 

9  o  )  .  (44) 

0 ;  Ybm 


eo  YbmR 


Next  the  charge  perturbation  6ois  obtained  by  carrying  out  the  momentum  space 
integration 

6p  =  -e  J <5f  d3p  ( 

to  find  the  fina^approximate  form  of  the  wave  equation 
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where 


(II  3  l2\.r  . 

\r  3 r  3r  r2/otz 


R.  iu2t>  /  gjBk\  f.u)A  /r-R  \ 

r  sf  l1  -  kr)  5Ez<r)  ~  *MRo>  -  ir  (it)  6Ez(r°} 


iou  yr-Rc 
o 


5{r-RQ+a)  5{r-RQ-a) 


“’pb^©!*2-  (r -R0V 


I,  wfJb\[  k2  SEz(Ro)  .  u?-2  6Ez(Ro)  tuo  (36Ez\  1 

\  kc  /  "  yT  SI2  W~  ft2  "  R  (4  \  3  r  L  Q 


os 


=  O)  -  iueQ  -  kBbc 


(47) 


(48) 


and  ®(x)  in  the  Heavyside  step  function 


®(x) 


x  >  0 
x<  0 


(49) 


Having  found  an  expression  for  5p  we  can  now  check  the  validity  of 

Equation  (30),  which  equivalently  neglects  the  6Jp  and  5J0  contributions.  The 

6J  contribution  is  of  order 
r 


(U(y2 « 1  1501 

and  the  6J0  contribution  is  of  order 
lV8 

~  «  1  -  (51) 
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THE  t-0  CASE 


We  shall  first  evaluate  the  i=0  case,  which  should  reproduce  the  well 
known  symmetric  space  charge  waves  of  the  hollow  beam.  The  wave  equation  reduces 
to. 


Ii_  r  !_  «  -  c 
r  3r  r  3r  otz  "  3 


5Ez(r)  -  <SEZ{R0) 


6(r  -  Rq  +  a)  <$(r  -  RQ  -  a) 


-  N  T  ^z(Ro)0|a2  -  (r  -  R0): 


(52) 


We  solve  this  equation  by  following  the  steps  from  Equation (88)  to  Equation  (39!  of 
Reference  8  to  find 


1  -  S 


where 


Kfc)  Kfe) (l  * 


-II1 


1  -  S  In 


ai 


R_  u>2 


s 

a  ws  ' 


) 


k2 


a? 


a 


S)  N  -  S 


•  (‘  - t) N  t ' 5 - s t N  t 

'  C  ’  O  C  CO 


(53) 

(54) 

(55) 

(56) 

(57) 
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After  Taylor  expansion  of  Equation  (53),  we  find 


(to  -  k3bc)2 


2v 

ZT 


.  3  .  1  2 

2£  2E 
1  +  e 


kc(kc  -  to3b)  In  — 


(58) 


where 


e 


“Bb\ 

wj 


and  v  is  Budker's  parameter  defined  by 

4irR  an.  e2 
_  o  b 


(59) 


(60) 


and  Nb  =  4TrRQanb  is  the  number  of  particles  per  unit  axial  length. 


When  e-K),  then  Equation  (58)  reduces  to 


(to  -  k3bc)2  =  h  (kV 
0  yb 


(61) 


g 

which  is  identical  to  the  result  obtained  by  Briggs  for  the  hollow  beam  and  solid 
beam  space  charge  waves.  In  obtaining  Equation  (61),  use  has  been  made  of  the 
approximation  to  =  <Bbc. 
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THE  1# 0  CASE 


We  solve  Equation  (47)  for  the  general  case  of  £.^0  by  following  the 
steps  from  Equation (88)  to  Equation  (99) of  Reference  8.  In  order  to  keep  the  size 
and  form  of  the  resulting  algebraic  eigenvalue  solution  short  and  simple,  we 
utilize  first  the  following  notation: 


c  _  Ro^b /,  ^b\ 

Si  -  rtl1  -  k r] 


.  %b(,  uSb\^o 

Sz  *  ‘  uf  V  '  kT  /  ~ 


Qz(Z2  -  1) 


.  ,  I,  “V 

^o  ^pbr  ~  kc 
n  'll 


The  eigenvalue  solution  becomes 


J4il2gf  +  2lS\  +  2^^)^+  (2igf  +  S,)  $x  1  -(§ jf*  } 

/  n  ,1  (R0\l 

x(l  +  Nj  +  N2)  +  (2tgf  +  Si)  (Xx  -  X3)]f  J  -  (24  +  Si)(X2  +  X„ )\^j 
+  s j(x2  -  x,)(^)4-  [24(gf  -  1)  +  S1(X1  +  x,)](^) 


+  22.NZ 


'BiA*  ffb 

y*z  \K 


\d  - 1  -  Si) (Si  -  so  +  r)  u  -  1  +  agf  +  so 


(Si  +  s2)  =  o 
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Equation  (77)  is  the  final  form  of  the  dispersion  relation  for  the 
modes  of  interest.  It  has  been  solved  numerically  for  the  growth  rate  =  In  f. 
and  the  Doppier-shifted  real  oscillation  frequency  nr  =  Re fi,  where  fi  =  w-kS^c  - 
for  a  broad  range  of  system  parameters,  y^,  <jjpb/ojc ,  wpt}R0/c,  a/RQ  and  R0/Rc> 

Typical  solutions  are  shown  in  Figures  3  for  yp  =  1.1  and  4  for  yb  =  3.  First 
of  all  we  note  that  for  azimuthal  mode  numbers  i =  0,  1,  and  2  the  solution  of  the 
dispersion  relation  of  Equation  (77)  gives  only  stable  waves.  Unstables  modes 
appear  for  l  &  3  as  shown  in  Figures  3  and  4.  Several  points  are  noteworthy  in 
Figures  3  and  4.  First,  and  contrary  to  all  previous  studies  of  the  diokotron 
instability,  we  observe  that  as  k-K)  the  diokotron  instability  disappears  and  is 
replaced  by  two  fast  waves.  For  low  y  beams.  Figure  3,  we  note  that  the  well 
known  f ilamentation  instability  occurs  at  small,  but  non-zero  wavenumber  k,  it  has  a 
spectrum  in  k  space  which  is  nearly  symmetric  about  the  k=0  point,  and  excludes  the 
k=0  point.  The  doppler  shifted  frequency  of  oscillation,  given  by 


is  constant,  to  lowest  order,  within  the  spectrum  of  instability  for  low  y  beams, 
y-1,  while  the  instability  growth  rate  decreases  away  from  the  maximum  occurring 
near  the  k-*±0  points.  Both  the  growth  rate  and  k-space  spectral  width  of  the 
instability  increase  with  the  azimuthal  mode  number  l.  The  real  frequency 
also  increases  with  i.  For  higher  y  beams  (Figure  4,  y^  =  3)  we  find  that  the 
lowest  order  symmetry  about  the  k=0  point,  that  was  observed  in  Figure  3  for 
yb  =  1.1,  new  disappears.  We  note  in  comparing  Figures  3  and  4  that  the  growth 
rate  and  the  doppier-shifted  real  frequency  decrease  when  the  y  of  the  beam 
increases. 


Figure  5  shows  the 
controlling  the  instability, 
we  note  that  both  the  growth 
applied  magnetic  field  while 
wavelengths. 


effect  of  varying  the  applied  magnetic  field  in 
For  the  positive  k  part  of  the  instability  spectrum 
rate  and  spectral  width  decrease  with  increasing 
there  also  occurs  a  shift  of  the  spectrum  to  longer 


2i 


NSWC  TR  81-385 


Figure  6  illustrates  the  dependence  of  the  growth  rate  on  the  beam 
energy  over  the  positive  k  spectrum  of  instability,  for  the  indicated 
parameters.  The  growth  rate  of  the  instability  is  monotonical ly  reduced  by 
increasing  the  beam  energy.  On  Figure  6  we  define  the  parameters  5  ,  and 
d?  for  convenience  in  further  analysis  and  discussion.  The  symbol  c0  stands  for 
the  long  wavelength  cutoff  limit  of  the  instability  positive  k  spectrum  and  the 
pair  of  symbols  d!7ijpb,  stand  for  the  value  and  associated  wave  number  of 
maximum  growth  rate.  Figure  7a  shows  the  variation  of  5  and  ^  with  yb  and 
Figure  7b  shows  the  variation  of  with  for  the  indicated  parameters. 

The  wave  numbers  5  and  increase  sharply  for  ^  1  and  quickly  saturate  to 
their  asymptotic  values  for  Yb^  1.5.  After  a  careful  examination  of 
Equations  (76)  and  (77)  we  find  analytically  that  the  asymptotic  value  of  r 
is  given  by 


k  c  u>.  . 

r  =  =  *6  -fit  i- 

co  w.  ^b 

pb  co 


(79) 


Figure  7b  shows  that 


rff  -  1/yJ  (80) 

in  agreement  with  the  basic  scaling  for  the  diokotron  mode  obtained  by  Uhm  and 
Siambis  in  Reference  1. 
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DISCUSSION  AND  CONCLUSIONS 


We  have  carried  out  an  extensive  analysis  of  the  filamentation  in¬ 
stability  of  hollow  beams  which  for  the  first  time  provides  a  detailed  analytical 
understanding  of  the  spectacular  experimental  results  obtained  by  Kyhl  and 
Webster  26  years  ago.  The  analytical  results  from  this  work  for  the  filamentation 
instability  when  compared  to  our  results  for  the  diokotron  instability  in 
Reference  1  show  that  the  two  instabilities  are  very  different  in  their  detailed 
properties  and  that  the  diokotron  instability  ( k=0 )  is  not  a  limiting  case  of 
the  filamentation  instability  when  k-O,  despite  the  fact  that  the  physical 
mechanisms  for  the  two  instabilities  are  of  a  similar  nature,  namely  the  azimuthal 
drift  of  the  electrons  and  the  hollowness  of  the  radial  density  profile.  This 
is  understandable  given  the  different  physical  situations  where  each  of  the  two 

instabilities  was  first  observed.  The  diokotron  instability  was  first  found  by 
2  5 

Buneman  ’  in  1944  as  an  explanation  for  magnetron  start  up.  In  France,  Doehler, 

5 

Warnecke  and  Mourier  supported  E-cross-B  flow  between  non-emitting  electrodes  and 
observed  the  diokotron  instability.  They  called  their  device  the  "diokotron" 
because  they  realized  the  effect  was  due  to  velocity  shear,  i.e.,  the  pursuit 
of  electron  layers  of  each  other  (ciuxeiv  *  pursue).  The  diokotron  instability, 
therefore,  found  birth  in  cross-field  devices,  most  prominent  of  which  is  the 
magnetron.  In  these  devices  the  electron  zero  order  motion  is  in  the  transverse 
(r,0)  plane,  where  the  diokotron  instability  also  takes  place.  The  axial  (z). 
direction  in  cross-field  devices  is  basically  an  ignorable  coordinate;  it  is 
often  short  and  wave  activity  along  it  is  suppressed  by  the  axial  boundary 
conditions  of  the  beam  and  the  confining  cavity,  hence  the  assumption  of  k=0 
for  the  diokotron  instability.  This  contrasts  very  sharply  with  the  physical 
situation  where  the  filamentation  instability  applies.  Namely  the  case  of  a 
hollow  beam  in  zero  order  axial  flow  having  only  a  higher  order  azimuthal  drift 
velocity  resulting  from  cross-field  forces.  In  this  case  wave  motion  along  the 
zero  order  beam  flow  in  the  z-direction  naturally  plays  a  significant  role;  hence 
k?0  and  the  filamentation  instability  has  a  spectrum  in  k-space,  in  fact  two 
spectral  bands,  one  for  positive  and  one  for  negative  k.  It  is  interesting  to 
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note  in  Figures  3,  4,  and  6  that  as  the  bean's  axially  directed  zero  order 
motion  decreases  to  zero  (>^-*1)  then  the  beau’s  azimuthal  drift  velocity, 
resulting  from  cross-field  forces,  becomes  more  important  and  the  fi lamentation 
instability  asymptotically  reaches  and  includes  in  a  limiting  sense  (as  k-*0) 
the  diokotron  instability. 


Fig.  6  Dependence  of  the  growth  rate  and  the  long  wavelength  cutoff  of  the  positive  k  spectrum  of  the 
instability  on  the  beam  energy  7^,  for  the  case  of  2  =  3,  a/R0  =  0.05,  R0/Rc  =  0.8,  Wpjj/cjc  =  0.3, 
cjpbR0/c  »  0.3.  (Note  that  the  latter  two  parameters  are  different  from  those  in  Figs.  3  -  5) 
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